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The Lieb-Oxford bound, a nontrivial inequality for the indirect part of the many-body Coulomb 
repulsion in an electronic system, plays an important role in the construction of approximations in 
density functional theory. Using the wavefunction for strictly-correlated electrons of a given density, 
we turn the search over wavefunctions appearing in the original bound into a more manageable 
search over electron densities. This allows us to challenge the bound in a systematic way. We find 
that a maximizing density for the bound, if it exists, must have compact support. We also find 
that, at least for particle numbers N < 60, a uniform density profile is not the most challenging for 
the bound. With our construction we improve the bound for N = 2 electrons that was originally 
found by Lieb and Oxford, we give a new lower bound to the constant appearing in the Lieb-Oxford 
inequality valid for any N, and we provide an improved upper bound for the low-density uniform 
electron gas indirect energy. 


I. INTRODUCTION 


A. Notation 


Lieb and Oxford (LO) [TJ |2] proved a nontrivial in¬ 
equality for the indirect part of the electron-electron in¬ 
teraction energy (total expectation of the interaction mi¬ 
nus the Hartree term) with respect to the LDA exchange 
functional. This inequality has been recently extended 
to include the gradient of the density [3]. The LO bound 
has played and continues to play a very important role 
in the construction of approximate exchange-correlation 
(xc) density functionals SHE]. While traditionally only 
the more general LO bound, valid for any number of par¬ 
ticles N (and corresponding to IV —> oo) has been taken 
into account in the construction of xc approximations, it 
has been shown very recently that the bound for N = 1 
and N = 2 is important in the context of metaGGA 
functionals mm, and can be imposed as an additional 
exact condition. 


The bound for N = 1 was first given in Ref. Q3J and 
proved rigorously in Ref. (2l For N = 2, Lieb and Ox¬ 
ford [2] could only provide a non optimal estimate of the 
constant appearing in the bound. In this work we de¬ 
velop a strategy to systematically challenge the original 
LO bound for a given number of electrons N. We use 
optimal trial wave functions for a given density, and we 
then vary the density in order to challenge the bound 
as much as possible. After showing that a density that 
maximally challenges the bound, if it exists, must have 
compact support, we follow the functional derivative of 
the bound to challenge it as much as possible without vio¬ 
lating Wrepresentability also for densities whose support 
is the whole space. As a first application of this proce¬ 
dure, we improve the lower bound for N = 2 given by 
Lieb and Oxford, see Eq. ( |60| below. Our construction 
also provides an improved lower bound for the constant 
appearing in the Lieb-Oxford inequality valid for any N, 
see Eqs. (63)-(64), and an improved upper bound for the 
indirect energy if the low-density uniform electron gas, 
see Eq. (|6l|). 


In electronic density functional theory (DFT) one is 
interested in finding the ground-state energy and density 
of iV-electron systems with Hamiltonian 

N 

H = T + V ee + U ext , U ext = ^ «( rj ). (1) 

i-1 

T and V ee are, respectively, the universal operators of the 
kinetic energy (in Hartree atomic units used throughout 
the paper), 


1 N 

*-2Xk (2) 

i =1 1 

and of the interaction (Coulomb repulsion) energy be¬ 
tween the N electrons, 


. 1 " 

^ 


1 - Si 


2 M 11* - : 

1,3 = 1 


( 3 ) 


The function u(r), in contrast, is a non-universal but 
arbitrary attractive external potential required to bind 
the repulsive electrons. Most of the formalism will be 
carried out for general spatial dimension D = 2 and 3, 
r £ R- 0 , focussing later on D = 3 only. 

In the following, T denotes a correctly normalized 
and antisymmetrized, but otherwise arbitrary Welectron 
wave function (thus not necessarily eigenstate of 0). 

= ^(rpri, ...,r N a N ), (4) 


where a n are spin variables. By pq,, we denote the par¬ 
ticle density associated with T, 

M r ) = n Y1 


J d D r 2 ...d D r N \'b(rai,r 2 a 2 , ...,r N <j N ) 


( 5 ) 
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B. Indirect Coulomb energy 


In terms of the local-density approximation (LDA) 


The electronic interaction energy in the quantum state 
T, defined as the expectation 

<*|fce|*> > 0, (6) 


excludes the infinite self energies of the point electrons, 
see the factor 1 — 5ij in Eq. If the electrons were a 
classical continuous distribution of negative charge with 
density p^(r), their interaction energy would be U[p^], 
with the Hartree functional 


U\p] = 


1 

2 


d D i 


d D r' 


P(r)p(r') 

|r —r'| 


> 0. 


(7) 


The most severe error introduced by this classical contin¬ 
uum approximation is a spurious finite self-interaction 
energy included for each electron. This is particularly 
evident in the case N = 1, since for any normalized one- 
electron wave function 'h, we have (^'|t4 e |^') = 0, while 
U\pis] > 0. The indirect interaction energy W\^>] is de¬ 
fined as 


W[V] = <tf|Vee|*} - U\p»\. (8) 

For wavefunctions that are ground states of an N- 
electron hamiltonian 0 (or good trial wavefunction for 
it) W'[\E r ] is normally negative. However, for a given den¬ 
sity p, it is possible to construct wavefunctions for 
which IT[d'] is positive or even infinity [2, 15]. We em¬ 
phasize that U\p) is a density functional, while IT[d/] is 
a functional in terms of the wave function ’I'. 


C. Lieb-Oxford bound 

The quantity IT['F] is limited by the Lieb-Oxford (LO) 
bound, 


C D J d D r p^{v) 1+1 / D < W[^}. 


(9) 


Cd > 0 is the unknown minimum possible number that 
makes this inequality true for all wave functions 4/ in 
D = 2 or 3 dimensions. So far, it is rigorously known 
that C 3 < 1.6358 [IB] and C 2 < 481.28 [T7] . and it has 
been argued [IS], on physical arguments, that the two 
bounds can be tightened to C 3 < 1.44 and C 2 < 1.96. 
The assumption behind these latter conjectured values is 
that the tightest possible bound is provided by the indi¬ 
rect energy of the uniform electron gas in the low-density 
limit, which, in turn, is commonly identified with the 
Wigner crystal total energy. This latter assumption has 
recently been proven wrong for the 3D case by Lewin and 
Lieb 0 ■ The study presented in this paper will also raise 
doubts on the first assumption that a uniform density is 
really the most challenging case for the LO bound, after 
a suitable optimal wave function for each given density 
has been defined (see Sec. IV). 


-A D J d D r P (v) l+1 / D = E^ DA [p] (10) 


to the D-dimensional exchange energy, with the exact 
constants A 3 = § (f ) 1 / 3 « 0.739, A 2 = |(f ) 1/2 « 0.798, 
Eq. ([9]) reads 

m < a D , (ii) 


where we have defined 


A[*] = 


IU[T] 

W^lP*V 


Ad = 


Cn 

Ad 


( 12 ) 


Considering all antisymmetric wave functions 4/ in D di¬ 
mensions, we may write 


A d = supAfd']. (13) 

m-.D 


The above rigorous upper bounds for Cd correspond to 
A 2 < 603, A 3 < 2.215. (14) 


Considering wave functions 4/ —> N with a given particle 
number IV, we define 

A d(N) = sup A[\l/]. (15) 

(<Z-.D)^yN 


Lieb and Oxford [2j have proven that A 3 (N) is monoton- 
ically increasing with its integer variable N, 

As(iV) < A 3 (V + 1), lim A 3 (1V) = A 3 . (16) 


They have also proven that A 3 (l) = 1.4786 (which was 
given originally by Gadre et al. M) and they have found 
a lower bound for A 3 (2), 

A 3 (2) > 1.67. (17) 

These bounds in D = 3 for N = 1 and N = 2 have 
been recently used to improve a certain class of exchange- 
correlation functionals 1111121- 

In this paper we develop a general strategy to find im¬ 
proved lower bounds for A d(N) by challenging the Lieb- 
Oxford bound, i.e, by evaluating A[4>] with particularly 
efficient trial wave functions T. Notice that this is dif¬ 
ferent from what is usually called tightening the bound, 
which means finding improved upper bounds to A d{N). 

A new lower bound for A_d(V) (or, generally, for Ad) is 
rigorously obtained each time we find a wavefunction that 
gives the highest value ever observed for A [4/] (for a given 
TV, or in general). Until very recently, it was believed that 
a lower bound for A 3 is given by A 3 > 1.444/A 3 = 1.955, 
corresponding to the total energy of the bcc Wigner crys¬ 
tal in the classical jellium model. However, in the jellium 
model, one can only identify the total energy with the in¬ 
direct energy if the electronic density is uniform, exactly 
equal to the one of the positive background. Only in 
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this case the electronic Hartree term will be exactly can¬ 
celed by the electron-background and the background- 
background contributions to the total energy. Lewin and 
Lieb [3] have shown that in the 3D case trying to make 
this cancellation happen by taking a superposition of all 
the possible Wigner lattices to have a uniform electronic 
density, introduces a shift that does not disappear in the 
thermodynamic limit. Thus, the value 1.955 does not 
correspond to the indirect energy of any wave f unct ion 
and is not a valid lower bound for A 3 . In Sec. IV we 


report a new lower bound for general TV, by considering 
an optimal trial wave function for TV = 60, and we also 
report an improved upper bound to the indirect energy 
of the low-density uniform gas. 


II. THE DENSITY FUNCTIONAL A[p] 


p(r) oc 

A[p] 

p(r) oc 

A[p] 

e -10(r-iy 

1.499 

e -50 (r-l Y 

1.262 

(1+r)- 4 

1.562 

e-^ 

1.689 

(1 + r)~° 

1.637 

e~ r 

1.699 05 

(1 + r)" 1 ’ 

1.662 

r e~ r 

1.698 66 

(1 + r)- 7 

1.674 


1.700 97 

(1 + r)- 10 

1.687 

r 1/3 e~ r 

1.700 95 

1 — r, r < 1 

1.638 

r -3 , re[Ri,R 2 \ 

1.145 

droplet 

1.498 

cos(r), r < ^ 

1.627 


Table I: Values A[p] for some simple spherical two-electron 
trial densities p(r) in three dimensions (TV = 2, D = 3), ob¬ 
tained numerically from Eqs. (A91-( All 1 of Appendix |A| In 
the last two rows we consider densities with compact support: 
“droplet” corresponds to the case of a sphere of uniform den¬ 
sity [24], and the density proportional to r - 3 m has been 
evaluated for Ri = 10 3 and R 2 = 10 5 . [Atomic units are used, 
where r is a dimensionless radial coordinate.! 


Considering only those wave functions —> p (in D di¬ 
mensions) that are associated with a given particle den¬ 
sity p = p( r), we define the density functional 

A[p] = maxA[T]. (18) 

4'->p 


Writing Ny = f d D rp^,( r) for the electron number in 
the state T, we then have 


A['k] < A.[p^r] < A d{N\$i) < Ad. 


(19) 


A. SCE interaction energy 


More explicitly, 

<M/|W ee |M/> - C7[p] 


A[p] = max 

'5>—>p 


£x LDA M 


K? e CE [p] - m 

E^ DA [p\ 


( 20 ) 


with the SCE interaction energy of Appendix [A] 

V e s e CE [p] = min(^|t4 e |^). (21) 

'f —>p 


The acronym SCE [191 - 121] stands for “strictly-correlated 
electrons” and defines a state |'I , sce[ / o] | 2 , which is a very 
accurate trial wave function (actually a distribuition) for 
the maximizing one in Eq. ( |20| ), being exact in ID [35] f° r 
any TV, and in any dimension for TV = 2 J23j. The SCE 
state is detailed in Appendix [A~[ In other words: Out of 
all antisymmetric wave functions T that are associated 
with a given density p, the one that provides (or is very 
close to) the strongest challenge to the Lieb-Oxford bound 
is the SCE state |^Esoe[ p]| 2 - Consequently, since I4| CE [p] 
can be e valuated rigorously for a wide class of densities, 
Eq. (A5) in Appendix [Aj we no longer need to consider 
different trial wave functions \k, but only different trial 
densities p instead, 


As a preliminary step, we have used simple analyti¬ 
cal trial spherical densities to evaluate A[p] for TV = 2, 
reporting the results in Table [I] We see that the lower 
bound ( fl7| ) is readily improved to 

A 3 (2) > 1.70097. (23) 


There is no need for considering scaled densities p^(r) = 
f D p(f r), with various values of £ > 0 , since A[p^] = A [p]. 
see Eq. (A14) in Appendix [A] It is interesting to notice 
that, once the most challenging wave function for each 
given p( r) is used, the densities that give the highest val¬ 
ues of A[p] are quite surprising. For example, a density 


proportional to e —50 (’ x ) , consisting of a thin spheri¬ 
cal shell, is similar to the one of the strongly-correlated 
limit of the Hooke’s atom series. Yet, it gives a value of 
A[p] which is much lower than the one obtained from the 
exponential density. Indeed, the strong-correlation limit 
of the Hooke’s series is known to give A[p] = 1.489 [S], 
again much less than what we obtain for exponential-like 
densities. The point is that previous works which ana¬ 
lyzed numerically the LO bound EM] focussed on physi¬ 
cal lramiltonians of the kind |l]), choosing v(r) that could 
be particularly challenging for the bound. In that con¬ 
text, exponential-like densities would correspond to the 
large nuclear-charge limit of the He isoelectronic series, 
which is a weakly correlated system. With our construc¬ 
tion, instead, we use the most challenging wave function 
for any given density, finding the unexpected trends of 
Table [T] We also see that the density of a uniform sphere 
(“droplet”) is not particularly challenging for the bound, 
a feature that will be further analyzed in Sec. |IV| for 
larger TV. 


B. Absence of a maximizing density without 
compact support 


X D = sup A[p], 

p:D 


Ap (TV) 


sup A[p]. (22) 

(, p:D)->N 


We now demonstrate that a function p(r) that max¬ 
imizes the functional A[p] for a finite TV cannot be a 
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physical density, unless it has compact support. The ar¬ 
gument is essentialy the same used by Lieb and Oxford 
[2] for N = 1 and N = 2. In terms of the SCE external 
potential of Appendix |A} 


where the indirect SCE interaction energy E e | CE [p] —C/[/)] 
is replaced with the functional E xc [p\ of the exchange- 
correlation energy, since the xc potential for IV-electron 
systems has the same asymptotic behavior as Eq. (29), 


MM _ ^e S e CE W 
vscvlp)(r) = 


and the Hartree potential 
I1M 5U[p] 

= JM 


= / d 


n r' 


p{r') 


(24) 


(25) 


Wc[p](r) = 


5E xc [p\ 


Sp(r) 


(r -> oo). (33) 


Quite interestingly, a density of exactly the same form of 
Eq. (30) for the 3D case, but only restricted in a finite 


we consider the Euler equation for maximizing A[p]. By 
writing p(r) = p( r) 2 to ensure p(r) > 0, and by varying 
p( r) we obtain 


vsce [p\ (r) - u H [p\ (r) Ke CE [p\ - U[p\ 


E EDA [p] 


E^ A \P} 2 


W W(r)?P(r) 


= PP{ r). 

If p( r) Y 0 everywhere, we obtain the Euler equation 

VsCE[p](r) -u H [p](r) 


(26, 


<5A[p] 

Sp(r) 


E™ A [p\ 

K! CE M “ U[p\ LDA 


e^W 


[p]( v)=p, (27) 


where p is the Lagrange multiplier ensuring fixed particle 
number N = J d D rp(r), and 


«* [p](r) = 


^ DA \p] 

Sp{ r) 


= -A 


D 


( 1 + 5>M 1/D 


(28) 

In this case, since u S CE[p](r) -A A-l and v H [p](r) -A A 
for r = |r| —>■ oo, we have asymptotically 


^SCE[p](r) - uh[p](i 


1 


(r -A oo). (29) 


Comparing this with Eq. (281, we see that a solution p{ r) 


of Eq. ( |27| must display the asymptotic behavior 

(r —> oo), 


, s h 

P( r ) -S’ rj 


(30) 


region of space (thus set to zero outside some region r € 
[i?i,i? 2 ]), has been considered by Perdew et. al. [TJ to 
study a general feature of GGA approximations related 
to the LO bound. Notice, however, that if we consider 
this kind of densities, p( r) a r~ 3 in r £ [i?i, A?. 2 ], even by 
choosing R\ and R 2 very large we get quite low values for 
A [p ], indicating that the asymptotic condition is anyway 
not enough to give a large A value, see Table]]] 

Even more generally, in a fictitious universe where the 
lectron-electron repulsion is multiplied by a factor a > 
1 , the density functional of their xc energy is given by 


E x 


M = J Q d/3{<^[p]|y ee |^[p]) - U\p\). (34) 


Here, out of all antisymmetric wave functions 4' that are 
associated with the same density p, \p] is the one that 
minimizes the expectaion (4'|T-t-/3V r e e|’I')! for any number 
(3 > 0. Since the corresponding a-dependent xc potential 
has the asymptotic behavior 


^xc,a [p] (l*) — 


<5-E xc ,a [/>] 


a 


6p(r) 


(r 


OO 


we conclude that even for the functional 


K[p] = 




M 


E^ A \p\ ’ 


(35) 


(36) 


the maximizing function p( r) must have compact sup¬ 
port. Notice that A a= i[p] = A[p] and lima^oo A a [p] = 

A[p]- 

If p(r) = 0 for |r| > r 0 , we see that, in principle, a 
maximizing density in Eq. (26) could exist. However, 


with our numerical investigation we have always found 
larger values of A for densities with unbounded support. 


with some constant k\. Such a function is evidently not 
normalizable, since with the D-dimensional volume ele¬ 
ment d D r = k 2 r D ~ 1 dr and a radius R > 0 finite but 
large enough, we have 


f d D rp{r) = I 
d| r |>i? Jr 


k\k 2 

dr - = 00. 


(31) 


We emphasize that this reasoning also applies to the 
modified functional 


AM = 


E xc [p] 

E^ A [pY 


(32) 


III. FOLLOWING THE FUNCTIONAL 
GRADIENT OF A[p] 

Although A \p\ has no maximizing density p without 
compact support, the functional gradient SA/5p tells us 
how to increase the value A[p] (or challenge the Lieb- 
Oxford bound) systematically. Starting from an N- 
electron density p = p(r) with a high value A[p], we 
consider a small density variation, 


P( r) 


p(r) + ea(r), J d 3 ra( r) = 0. 


(37) 
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Provided that eer(r) is truly “smallwhich precisely 
means that 


/ d 3 r <r(r) 2 = 1 


(38) 


and |e| <C 1, we have 

A[p + ea] - A[p\ « e J d 3 r G[p](r) c?(r), (39) 

with the gradient G[p](r) = 5A[p\/8p(r) given by 
Eq. ( |27| . Although fd 3 rcr(r) = 0, the right-hand 
side of Eq. (39) can nevertheless be > 0, provided that 


G([p];r), as a function of r, is different from a constant, 
G[p](r) ^ const. 


A. Formal optimization of the increment 

Formally, maximizing the integral f d 3 r G[/?](r)cr(r) 
with respect to er(r) subject to the two constraints 
f d 3 rcr(r) = 0 and f d 3 rcr(r) 2 = 1, 

<Rr){ / d3rG W( r ) CT ( r ) J d3r,(J ( r ) 

-p 2 Jd 3 ra(r) 2 ^ = 0, (40) 

yields the Euler Equation G[p](r) — /xi — 2p 2 a(r) = 0, 
with the solution 


(r) = 


G[p\(r)- p 1 
2/i2 


(41) 


The first Lagrange multiplier pi is fixed by the normal¬ 
ization constraint f d 3 rao(r) = 0. The second one p 2 is 
absorbed in the small parameter e, guaranteeing the va¬ 
lidity of the approximation ( [39] ). [Independently, small¬ 
ness of ea( r) is necessary (but not sufficient) for the re¬ 
sulting density to be non-negative, p( r) + ecr(r) > 0 for 
all r.l 


For a A’-electron (hnite) density p , Eqs. (27)-(29) im¬ 
ply the large-?' behavior (r —> oo) 


G[p](r) 


1 


rl 


\ E * DA [p}\ Lr 


A[p\A D (l + ^ p(r) 1/D 


(42) 

Necessarily, er 0 ( r ) —> 0 for r —► oo, implying pi = 0 in 
Eq. plj ). Consequently, due to the term 1/r in Eq. (42), 
Jd 3 r<ro(r) cannot be zero (or even finite). In other 
words, p( r) + ecro(r), with e / 0, must, again, yield a 
density with compact support. In the following, we give 
an analytical example for the case of a density with com¬ 
pact support. 


B. Analytical example for densities with compact 
support 


As an example, we evaluate Eq. (41) for the spherical 
2-electron density [233 


P ( r ) = 


Po ( r < R), 
0 (r>R) 


Po ~2 ttR 3 ' (43) 


This density corresponds to a uniformly charged sphere 
(“droplet”) with radius R and total charge 2, 


U\ P \ = 


12 

5R’ 


vn[p]{r) = 


3 R 2 - r 2 
R 3 


(r < R). (44) 


The exact exchange energy E x [p\ is given by —E x [p\ = 

\ u \p\ = IT ? whiie m 


-E. 


LDA 


\P] = 


1.1545 
R ’ 


(45) 


^ LDA WW = -(^) 1/3 ^ ( r<R ). (46) 

From Ref. [24j . we have A \p] = 1.498 and the SCE co¬ 
motion function (see Appendix [A] is 


r 3 \ 1/3 


f(r) = R{l R ,j 

The resulting SCE external potential is given by 

» sc "H(r)=, SCE M(0)- f du 

Jo 


(47) 


f(uW 


= « SOE W(0)-ijf' 


~/R 


dx 


o [a; + (1 — x 3 ) 1 / 3 ] 2 


Eventually, Eq. (41) reads 

—?; SCE [p](r) + UH[p](r) - jh 


cr(r) = 


2jh 


(48) 


(49) 


where /A = p\ — A[p]?; EDA (note that f EDA does 
not depend on r in the present example) and 2fx 2 = 
—E x DA [p]2p 2 > 0. The constant u SCE [p](0) can also 
be absorbed by the multiplier which is fixed by the 
conditon J R dr 47rr 2 cr(r) = 0. Then, we have 

1 f r 2 r r/R dx 1 

a{r) = WiRi&~^ ll + L [x + (i-x 3 y / 3 ] 2 }' 

(5°) 

A simple but accurate approximation to this function (for 
R = 1) is 

1 


TapprW - 2 - 2 L 


0.4r 3 -1.85r 2 +r+0.16 r< 1.(51) 
We therefore consider the densities (for r < 1) 


Pa(r) = po+a 


0.4r 3 -1.85r 2 +r+0.16 


(a > 0) (52) 


to obtain the values A[p 0 2 ] = 1.521, A[p 0 5 ] = 1.551, 
Afpr.o] = 1.590, A[ Pl . 5 ] = 1.611, A[ Pl . 6 ] = 1.612. [For 
a > 1.6, the density p a (r) becomes negative.] 
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C. Compromise for A-representability for densities 
with unbounded support 


Since we observe higher values of A [p] for densities for 
which p(r) ^ 0 everywhere, we consider here a compro¬ 
mise to follow the gradient of A[p] without violating N- 
representability. We perturb the density with some func¬ 
tion cr(r) in Eq. (37) that depends on a certain number of 
parameters and satisfies Eqs. (37)-(38), keeping the per¬ 
turbed density A-representable with suitable constraints. 
One can then choose the parameter values in order to 
maximize the overlap with the gradient [the right-hand- 
side of Eq. (39)]. 

As an example, we start from the 3D exponential two- 
electron density (D = 3, N = 2) 


= g, (53) 

which already gives the high value A[p] = 1.69905 (see 

Table [I]). We choose for cr(r) the parametrized form 

CTa(r) = \/~( 1 ^y) e_ar ( 54 ) 

which obeys the conditions f d 3 rer a (r) = 0 and 

J d 3 rcr a (r) 2 = 1 for all values of the parameter a > 0, 
so that the function p a ,e( r ) = p(r) + e<r a (r) is always 
correctly normalized. In addition, A-representability re¬ 
quires that p a ,e(r) > 0 for all r > 0. For any value 
of a, this is fulfilled for e m i n (a) < e < e max (a), where 
£min(a) < 0 and e max (a) > 0 are given by 


emin(&) 


3( a ~1) (3— 4a)/a 

4aV3a 3 n 

1 


(0 < a < |), 

«>!) 


€ 

A [p + ea a2 \ 

A[p] + el(a 2 ) 

0 

1.699 052 

1.699 052 

-0.01 

1.700 487 

1.701 354 

-0.015 

1.700 833 

1.702 505 

-0.02 

1.700 868 

1.703 656 

-0.0207 

1.700 843 

1.703 817 


Table II: Exact values A[p + ea a2 ] for various values of e < 0, 
compared with the first-order expansion. 


47 t = —0.2302 at 02 ~ 2.49. While e(ai) = 0, we have 
£( 02 ) = —0.0207, and Eq. ( |58| ) for a = 02 gives 

A[p + e(a 2 )a a2 \ - A[p\ « 0.004 765. (59) 


In Table|H|we report the values of A[p+ecr a2 ] as a function 
of e and compare them with the ones from the first-order 
expansion. As predicted, we see that A[p] increases for 
small e. However, the first-order expansion breaks down 
before e(a 2 ), so that the maximum value of A[ p] t hat we 
obtain is less than the one predicted by Eq. (59). The 
improvement in this case is very small, but we suspect 
that this is due to the fact that for N = 2 the exact A 3 (2) 
is very close to 1.701, so that we are really hitting the 
boundary. In fact, in the previous example of Sec. |IIIB| 
we have seen that when we start from a much less optimal 
density the improvement in A[p] with our procedure is 
much larger. 

We have also repeated the procedure using as a starting 
density the one corresponding to e = — 0.02 in Table |Tl} 
but we could only slightly improve the result obtaining 
A[p] = 1.701052, which is, so far, our best value, 


A a (2) > 1.701052. (60) 


^max(^) 


3(a—1) ^—(3—4 a)/a 
4a\/3a 3 7r 


(0 < a < 1 ), 
(a > 1 ). 


(56) 


IV. IS A UNIFORM DENSITY THE MOST 
CHALLENGING FOR THE LIEB-OXFORD 
BOUND? 


Evaluating numerically the functional gradient 
G[p](r) = (5A[p]/<5p(r) of Eq. (27) for the density 
p(r) = p(r) of Eq. (53), we consider, as a function of a, 
the overlap integral 

/»oo 

1(a) = / dr( 47 rr 2 ) G[p](r) cr a (r). (57) 

Jo 

For any value of a > 0, the maximum possible value in 
Eq. (391 is approximately (if the first-order expansion 
holds) 


A[p + e(a)a a ] - A[p] « e(a)/(a), (58) 

where e(a) = e max (a) > 0 for 1(a) > 0 and e(a) = 
Emin (a) < 0 for 1(a) < 0. 

Numerically, 1(a) > 0 for 0 < a < 1 and 1(a) < 0 
for a > 1, with a strong maximum I(a{) ~ 2.9267 • 47 r 
at a\ « 0.079 and a weak minimum I(a 2 ) ~ —0.01479 • 


In Ref. |T5] it has been argued that the tightest bound 
should correspond to the case of the uniform electron gas 
at extremely low density (equivalent to the SCE limit 
for a uniform density). This suggestion was made by 
considering electronic hamiltonians of the form 0 with 
particularly challenging w(r), keeping in mind that the 
bound increases [2 with the number of electrons N, see 
Eq. @. 

With our formalism, we directly consider the most 
challenging wavefunction (or one which is very close to 
it, thus providing anyway a lower bound for A[p]) for 
each given density, and we can thus question whether a 
uniform density profile is really the most challenging for 
the bound. Already by putt ing together existing data, 
we can compare, in Table III the values of A[p] obtained 
from the (sphericalized) atomic densities of Li, Be, C, 
B, and Ne [21] , with the ones obtained from spheres of 
uniform density (“droplets”) [M] : we clearly see that the 
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TV 

A[p] atomic 

A[p] droplet 

3 

1.713 

1.550 

4 

1.731 

1.603 

5 

1.747 

1.627 

6 

1.767 

1.657 

10 

1.816 

1.708 


Table III: For different values of TV, we compare the A[p] ob¬ 
tained from atomic densities (values from [21]) with the ones 
obtained from spherical droplets of uniform density (values 
from [24]). 



Figure 1: Values of A[p] for the fixed density profile p(r) oc 
r 1 / 2 e -r compared with those for spheres of uniform density 
(“droplets”) as a function of the particle number TV. The val¬ 
ues for p(r) oc r 1 ' 2 e _r are significantly higher than those for 
uniform densities. The size extrapolation for uniform densi¬ 
ties is also shown, where the fitting parameters are a = 1.918, 
b = —0.3253, c = —0.2791. 


finding a = 1.918, b = —0.3253, c = —0.2791. The fitting 
function is also shown in Fig. |T] The value of the fitted 
parameter a gives our TV —► oo extrapolation for A[p] in 
the uniform electron gas, A uni f[TV —>• oo] = 1.918. This 
value can be compared with the one obtained by taking 
the r s — > oo limit of popular LDA parametrizations: for 
example, the PW92 |2B] parametrization yields 1.947 at 
zero polarization and 1.977 for the fully polarized case, 
while the VWN [ST) at zero polarizations gives 1.9043. 

After Lewin and Lieb [3] showed that the value 1.955 = 
1.4442/A3 does not correspond to an indirect energy, 
our value A[p] = 1.91175 for TV = 60 and spherically- 
symmetric density profile p(r) oc r v ! 2 e~ r is the highest 
value of AsfT] ever observed, setting a new lower bound 
for A 3 (TV) for any TV, so that, rigorously 

1.91175 < A 3 < 2.215, (63) 

or, in terms of the constant C 3 in Eq. ([9]) 

1.4119 <C 3 < 1.6358. (64) 


V. CONCLUSIONS AND PERSPECTIVES 

We have developed a method to maximally challenge 
the Lieb-Oxford bound, using optimal (or nearly opti¬ 
mal) trial wave functions that can be constructed from 
a given density. This allows us to rewrite the most chal¬ 
lenging bound for a given number of particles directly as 
a density functional. As a first application of the method, 


atomic densities yield significantly higher values of A [p]. 
as already observed for TV = 2 in Table [T] 

We have also performed calculations with the fixed 
spherical density profile p(r ) oc r 1 ^ 2 e~ r , which was par¬ 
ticularly promising for TV = 2 (see Table [T]), for particle 
numbers TV < 60, and compared the values with the ones 
for spheres of uniform density, extending the calculations 
of Ref. [211 up to TV = 60. The results are reported in 
Fig-E where we clearly see that the uniform droplets give 
values significantly lower for A[p]. This suggests that a 
similar behavior may arise in the limit TV —> 00 : a density 
with particular modulations might challenge the bound 
more than the uniform one. 

Our new value for the uniform sphere at TV = 60, A = 
1.818, sets an improved upper bound [3], equal to —1.343, 
for the low-density uniform electron gas indirect energy 
per particle w, which then must be between 

- 1.45 < w < -1.343, (61) 

where the lower bound —1.45 has been proven in (25] . 
We have also performed a size extrapolation of our A[p] 
for the droplets of uniform density of Fig. [l] by fitting 
our data to a lquid-drop model expansion 

Auni f [TV] =a + bN~^ 3 +cN~ 2 / 3 , (62) 


• we improved - see Eq. ( |60[ ) - the constant in the 
LO bound for TV = 2, which provides a constraint 
to develop new metaGGA functionals mi; 


• we have given an improved lower bound for the con¬ 
stant appearing in the LO inequality valid for all 
particle numbers TV, see Eqs. (63)-(64); 


• we have obtained an improved upper bound for the 
indirect energy per particle of the low-density uni¬ 
form electron gas, see Eq. (611. 


In future works we will analyze systematically the bound 
for larger particle numbers TV, trying to give improved 
lower bounds for A_d(TV) and for Xp. 

More generally, from this study we have learned that it 
is quite difficult to predict which densities will maximally 
challenge the bound (see for example Table [ij the trends 
reported there seem totally unpredictable). For sure, we 
observe that, for finite TV, a uniform density is not the one 
that challenges the bound the most, suggesting that the 
indirect energy of the uniform gas at low-density may 
not provide the tightest bound, contrary to what was 
previously suggested. 
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As usual, the functional derivative is determined up to 
a constant, which for finite systems we fix by requiring 
that the potential vanishes at infinity. 

For example, a spherical two-electron density p(r) = 
p(r) in D dimensional space has the co-motion function 

Pam 

f 2 (s) = -/( S )n- (A8) 

M 

In terms of the invertible function 

N e (s)= f d D rp(r), (A9) 

J|r|<s 


Appendix A: Strictly correlated electrons 


The minimizing antisymmetric wave function = 

d’KsIjo] hr the definition of the density functional of the 

non-interacting kinetic energy [28| . 

T s [p\ = min(*|T|*>, (Al) 

W —>p 

is usually a Slater determinant of Kohn-Sham orbitals. 
In contrast, the minimizing one in 

V e S e CE \p] = min (* ivy*) (A2) 


is (or it is very close to) a state ’P = TsceIp] with strictly 
correlated eletrons (SCE). v Psce[p] is not a regular wave 
function but a Dirac-type distribution. Its position rep¬ 
resentation is singular, 

4' SCE [p](ri<Ti,...,r A rcr Ar ) = 0 for (r 1; ..., r N ) £ fl 0 [p]- 

(A3) 

Here, Hq [p] is a D-dimensional subspace of the (N x D)- 
dimensional configuration space of the IV-electron sys¬ 
tem, given by 


^0 [p\ = {(s,f 2 (s),...,fiv(s)) 


s G S p C 




}> (A4) 


where S p = { re R 13 | p( r) ^ 0} is the spatial region of 
non-zero density. fi(s) = s, f 2 (s), ..., fjv(s) are co-motion 
functions: In an SCE state, a configuration (r 1; is 

observable only when its N positions obey the relations, 
r„ = f n (ri), n = 2Then, the distance between 
electrons i and j is |fj(s) — f, (s)|, fixed by the position 
iq = s of electron 1. Consequently, we have 


ySCE 

•'ee 


\p\ 



N— 1 N 


E E 

i—l j=i +1 


i 

i f i(s) -fj(s)r 


(A5) 


This is truly a density functional, since the co-motion 
functions are fixed by the density, f n (s) = f„[p](s). For a 
large class of densities p, the functions f„[p](s) and thus 
the functional I4e° E [p] can he evaluated rigorously na 
I2TK231 . Its functional derivative turns out to be [29} f30j 

= t>scE[p](r), (A6) 

with the SCE external potential nscE[p](r), fixed by 


Vn S CE[p](r) 


N 


-E 


r - f ra (r) 
l r -f n ( r )| 3 ' 


(A7) 


the radial co-motion function is given by 

f(s) = N~ 1 (2-N e (s)). (A10) 


Eq. 

Eq. 


(A51 


(A8| implies that jr 


f 2 (r)| = r + f(r ) and, due to 


ySCE 

^ee 


\p\ 


1 

2 


d D r 


P(r) 

r + f(r)' 


(All) 


Due to Eq. (A7), the SCE external potential, with 
ngcE —> 0 for r —> 00 , is 


^SCE[p](r) 


ds 


+ f(sW 


(A12) 


For any A-electron density p(r) with co-motion func¬ 
tions f n (r) (n = 2,..., A), we may consider the continu¬ 
ous series of scaled 7V-electron densities p^(r) = ^ D p(£r ), 
with ^ > 0 and f d D rp^(r) = f d D r p(r) = N. The 

co-motion functions (r) of pj(r) are given by 


f ™ ?) ( r ) = ^fn(?r)- 


(A13) 


Therefore, the functional (A5) has the simple scaling 
property 


K s e CE [pd = nv£y*[p}. 


SCEr 


(A14) 


We should remark that the SCE wave function as a 
minimizer for the electron-electron interaction energy has 
been first conjectured on physical grounds [13 EH- In re¬ 
cent years, it was reco gnized that the problem posed by 
the minimization (A21 is equivalent to an optimal trans¬ 


port problem with Coulomb cost E3EU- Since then, 
the optimal transport community has produced several 
rigorous results. In particular, the SCE state has been 
proven to be the true minimizer for any N in ID [22] 
and in any dimension for N = 2 [231 . For more general 
cases, it has been shown that the minimizer might not 
be of the SCE form [32] . Even in that case, however, 
SCE-like solutions seem to be able to go very close to the 
true minimum 1331. 
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